m  •  <  .'uo'0«b  jozti 


2-  jj  •  1  ij •  ;i /  AM  iWA  rro-rsms.  MiO-W- 


REPORT  DOCUMENTATION  PAGE 


AFRL-SR-BL-TR-98- 


t  iwwwt  mum  w  i 


1.  AGENCY  USA  OMLT  «***• 


Mtttx*  Of<kra**MM 

••  tss  rr  sbsst  sztT 

7/15/98 


^  jjjj  i  AHp  VllTtni 

"High  Frequency  Electromagnetic  Propagation/Scattering 
Codes" 


Fi^l^fethniSI  TifSr'r^/97-5/14/98 


1  FUMOUia  NUMttHS 

pf%2o-  97- C- 00*53. 


«.  aUTNOMS) 

Vladimir  Rokhlin 


rT.T.Tii^fi4i 


t  rrwntMWB  *w 

Fast  Mathematical  Algorithms  &  Hardware  Corporation 
1020  Sherman  Avenue 
Hamden, CT  06514-1337 


a.  ACRPOIIMIfM  owgani&atnm 
HPOR?  MlllMIR 


I.  jPONSOWNGyMUMOTOMwe  AGENCY  MAMC(SI  ANO  AOOUSHE. 

USAF/AFMC  . 

Air  Force  Office  of  Scientific  Research/A^' 
110  Duncan  Avenue,  Room  B115 
Bolling  AFB,  DC  20332-8050 

II.  jijr*U*ttNTA»Y  MOTES 


10.  SPONSORING  /  MOM 

agency  wort  numsu 


ia*.  aanuaunoN  /  a*ar  ahiuty  >tat 


12b.  OtSTRIBUTlOH  COOC 


Distribution  is  unlimited;  approved  for  public  release 


^^BJ^!7i5«iSnM5nrffr:mrr*y|lPrpnrt  «i°vQlr>porl  1  indc‘T>  ^ttr  Contract" 

l^enever  bsnd-limited  signals  are  measured  or  generated,  different  ctistrituticns  of  a  fixed  nuiter  of 
receivers  and  traisdLcers  lmd  to  very  different  resolutions.  A  closely  related  set  of  issues  is  encantered 
in  tie  numerical  solution  of  scattering  problems;  given  a  scatterer,  one  would  like  to  find  nodes  cn  its  air- 
face  leading  to  most  efficient  discretizations.  During  Fhase  I  of  the  Si  Ik  project  F®629-97-O0Q52,  we  dis¬ 
covered  (somewhat  werendipitoulsy)  that  whenever  band-limited  signals  are  to  be  discretized,  measured,  or  ^ 
^TPrgfpri  the  ccretructicn  of  optimal  (in  a  very  strong  sense)  dcnfigjratirns  of  nodes  is  a  tractable  problem. 
WiEn  the  rubs  are  to  located  cn  a  line  or  on  a  disk  in  FT,  the  solution  is  a  fairly  stiairjitfarward  ccnsecpen 

of  classical  reseats  dbtaired  by  Slepian  and  his  collaborators  mare  tten  3D  years  ago.  The  reqiired  arrerioal 

tmi«;  are  grit*  efficient;  the  resulting  discretizations  are  a  dramatic  inprovement  over  the  ones  currently 

rt~im  of  optimal  discretizations  of  more  ccnplicated  regicnssrecpnies  additional  mathematical^ 
apparatus;  such  spparatus  tas  been  largely  (thou^i  by  no  means  conpletely)  desiged.  liTfartunately,  the  runer 
cal  tools  we  have  ccnstiueted  are  gate  inefficient  in  this  environment,  which  limits  the  size  of  regions  on 
vWch  we  are  currently  ahla  to  treat  effectively.  We  have  applied  far  the  Fhase  n  of  theSTIR  project  to 
further  develop  such  analytical  and  numerical  tools,  and  to  apply  them  both  in  the  engineering  and  nunerical 

environments. 


ITT  WTTTTTT  wf^n 


5 ttK,  i , 


“High  Frequency  Electromagnetic 
Propagation/Scattering  Codes” 


Final  Technical  Report 

July  15, 1998 


Small  Business  Technology  Transfer  (STTR)  Program  (DOD) 

DCMC  Hartford 

Topic  No:  AF97T001 

Issued  by  USAF,AFMC 
Air  Force  Office  of  Scientific  Research  under 

Contract  No:  F49620-97-C-0052 


F.M.A.&  H.  Corporation  P.I.:  Paolo  Barbano,  Res.Sci. 

1020  Sherman  Avenue  (203)248-8212 

Hamden,  CT  06514  This  report  covers  period: 

Eff.  Date  of  Contract:  8/15/97  8/15/97  -  5/14/98 


“The  views  and  conclusions  contained  in  this  document  are  those  of  the 
authors  and  should  not  be  interpreted  as  representing  the  official  policies, 
either  expressed  or  implied,  of  the  Small  Business  Technology  Transfer 
(STTR)  Program  or  the  U.S.  Government.” 


“Distribution  is  unlimited;  approved  for  public  release.” 


DTIC  QUALITY  mCPBCTED  1 


High-Frequency  Electromagnetic  Propagation/ Scattering 

Codes 


1  Introduction 

Whenever  band-limited  signals  are  measured  or  generated,  the  locations  of  receivers 
or  transducers  have  to  be  selected;  it  is  well-known  that  different  distributions  lead  to 
very  different  resolutions  given  a  fixed  number  of  receivers  or  transducers.  A  closely 
related  set  of  issues  is  encountered  in  the  numerical  solution  of  scattering  problems: 
given  a  scatterer,  one  would  like  to  find  nodes  on  its  surface  leading  to  most  efficient 
discretizations. 

During  Phase  I  of  the  STTR  project  F49620-97-C-0052,  we  discovered  (somewhat 
serendipitously)  that  whenever  band-limited  signals  are  to  be  discretized,  measured,  or 
generated,  the  construction  of  optimal  (in  a  very  strong  sense)  configurations  of  nodes 
is  a  tractable  problem.  When  the  nodes  are  to  be  located  on  a  line  or  on  a  disk  in 
R.2,  the  solution  is  a  fairly  straightforward  consequence  of  classical  results  obtained  by 
Slepian  and  his  collaborators  more  than  30  years  ago.  We  have  constructed  the  necessary 
numerical  tools,  which  are  quite  efficient;  the  resulting  discretizations  are  a  dramatic 
improvement  over  the  ones  currently  employed. 

Construction  of  optimal  configurations  of  nodes  on  more  complicated  regions  requires 
additional  mathematical  apparatus;  such  apparatus  has  been  largely  (but  by  no  means 
completely)  designed.  Unfortunately,  the  numerical  tools  we  have  constructed  are  quite 
inefficient  in  this  environment,  which  limits  the  size  of  regions  on  which  we  are  currently 
able  to  design  optimal  configurations  by  approximately  three  wavelength.  Needless  to 
say  that  this  state  of  affairs  is  not  satisfactory;  we  are  currently  working  on  improved  (in 
terms  of  CPU  time  requirements)  algorithms  for  the  construction  of  such  configurations. 

At  the  present  time,  we  are  finalizing  the  first  of  a  series  of  papers  describing  this 
work  (this  paper  describes  the  one-dimensional  version  of  the  theory).  We  have  filed  a 
patent  disclosure  with  Yale  University  (of  which  we  are  attaching  a  copy),  and  applied 


for  the  Phase  II  of  this  STTR  project.  The  present  report  describes  the  state  of  affairs 
as  of  the  conclusion  of  Phase  I,  with  an  emphasis  of  the  analytical  and  numerical  (as 
opposed  to  the  engineering)  aspects  of  the  situation.  It  should  be  kept  in  mind  that  the 
results  discussed  here  are  very  recent,  and  our  view  of  the  field  is  changing  constantly. 


2  Background 

When  measurements  are  performed,  it  often  happens  that  the  signal  to  be  measured  is 
well  approximated  by  linear  combinations  of  oscillatory  exponentials,  i.e.  functions  of 
the  form 


;=i 


e''X)'x 


(1) 


in  one  dimension,  of  the  form 


X>. 

i=i 


(2) 


in  two  dimensions,  and  of  the  form 

n 

Y^otj  •  e‘-(Vx+*VS'+t/j  -2) 
j= 1 


(3) 


in  three  dimensions.  In  most  cases,  the  signal  is  band-limited,  i.e.  there  exist  such  real 
positive  a  that  all  1  <  j  <  n, 


I  Aj  |<  a 

in  one  dimension, 

+  fx ]  <  a 2 

in  two  dimensions,  and 

+  /Xj  +  v2  <  a2, 


(4) 


(5) 


(6) 
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in  three  dimensions. 

As  is  well-known,  most  measurements  of  electromagnetic  and  acoustic  data  (espe¬ 
cially  at  reasonably  high  frequencies)  are  of  this  form.  Examples  of  such  situations 
include  geophone  and  hydrophone  strings  in  geophysics,  phased  array  antennae  in  radar 
systems,  multiple  transceivers  in  ultrasound  imaging,  and  a  number  of  other  applications 
astrophysics,  medical  imaging,  non-destructive  testing,  etc. 

In  this  proposal,  we  describe  a  procedure  for  determining  the  optimal  distribution  of 
sources  and  receivers  that  maximizes  accuracy  and  resolution  in  measuring  band-limited 
data  given  a  fixed  number  of  receivers.  Alternatively,  the  procedure  can  be  used  to 
determine  the  optimal  distribution  of  receivers  that  will  minimize  their  number  given 
specified  accuracy  and  resolution.  We  illustrate  the  procedure  with  several  examples 
of  optimal  distributions  of  receivers  on  an  interval  (one- dimensional  case),  on  the  disk 
(two-dimensional  case),  and  on  the  triangle  (two-dimensional  case). 

Remark  2.1  Our  procedure  permits  the  design  of  optimal  receiver  distributions  in  ar¬ 
bitrary  regions  in  one,  two,  and  three  dimensions.  However,  when  the  region  in  question 
is  an  interval,  a  disk,  or  a  sphere,  there  exists  an  explicit  formula  for  the  locations  of 
the  receivers,  in  terms  of  zeroes  of  the  Generalized  Prolate  Spheroidal  Wave  Functions. 
Somewhat  more  complicated,  but  still  nearly  analytical,  schemes  exist  when  the  region 
is  the  triangle  (in  two  dimensions),  a  tetrahedron  (in  three  dimensions),  and  in  several 
other  cases.  In  order  to  avoid  describing  the  numerical  apparatus  required  in  the  general 
case,  we  restrict  our  examples  to  the  cases  where  the  problem  can  be  solved  explicitly 
(as  opposed  to  an  involved  numerical  procedure). 

The  analytical  apparatus  we  use  is  closely  related  both  to  classical  Gaussian  quadra¬ 
tures  and  to  the  analysis  of  band- limited  functions  initiated  in  the  classical  works  [l]-[4]. 
It  results  in  optimal  (in  a  very  strong  sense)  choice  for  both  source  and  receiver  loca¬ 
tions.  While  the  construction  is  quite  straightforward  in  one  dimension  (in  this  case,  it 
is  a  remarkably  simple  corollary  to  [1],  [2]),  in  higher  dimensions  it  is  complicated  by 


the  fact  that  regions  can  be  of  various  shapes.  In  practical  measurements,  the  cases  of 
the  string,  of  the  disk,  of  the  square,  and  of  the  triangle  appear  to  be  most  frequently 
encountered. 

3  Outline  of  Theory  In  One  Dimension 

Given  a  real  a  >  0,  we  will  be  considering  the  operator  Pa  :  L2[— 1,1]  — >  L2[— 1, 1]  defined 
by  the  formula 

PM(x)  =  £  ef— •  V(t)dt,  (7) 

for  any  tp  €  L2[— 1,1];  the  image  of  L2  under  Pa  will  be  denoted  by  Xa.  As  observed 
in  [1],  all  eigenvalues  A0,  Al5  A2,  •  •  •  of  the  operator  Pa  are  distinct,  and  the  corresponding 
eigenvectors  are  the  appropriately  chosen  Prolate  Spheroidal  Wave  functions. 

The  following  theorem  is  a  combination  of  several  Lemmas  from  [1]. 

Theorem  3.1  For  any  positive  real  a,  the  eigenfunctions  of  the  operator 

Pa  constitute  an  orthonormal  basis  in  Xa.  The  function  fjj  has  exactly  i  zeroes  on  the 
interval  [—1, 1],  and  for  any  positive  integer  n,  the  functions  ip0,  ^i?  •  •  • ,  f’n  constitute  a 
Chebychev  system  on  the  interval  [—1,1]. 

Assuming  that  A0,Ai,A2,***  are  ordered  so  that  \  A0  |>|  Ai  |>|  A2  |>  •  •  •,  the  even- 
numbered  eigenvalues  are  purely  real,  and  the  odd-numbered  ones  are  purely  imaginary. 
The  even-numbered  eigenfunctions  are  even,  the  odd-numbered  ones  are  odd,  and  the 
j  —  th  eigenfunction  has  exactly  j  roots  on  the  interval  [—1,1]. 

Finally,  for  all  sufficiently  large  i, 

I  ^  |<  e->.  (8) 

Remark  3.1  The  operator  (7)  and  some  of  its  generalizations  have  been  studied  in  [1]- 
[5];  later,  a  systematic  investigation  of  generalizations  of  (7)  was  undertaken  in  [7]-[9]. 
The  investigations  reported  in  [l]-[5]  appear  to  have  been  motivated  almost  entirely  by 


considerations  of  Electrical  Engineering;  the  work  [7]-[9]  seems  to  have  been  motivated 
predominantly  by  considerations  of  pure  mathematics.  Remarkably,  implications  of  this 
mathematical  apparatus  in  numerical  analysis  appear  to  have  been  overlooked. 

In  the  interests  of  brevity,  we  will  be  referring  to  the  eigenfunctions  tpn  of  the  oper¬ 
ator  (7)  as  Slepian  functions,  and  to  their  roots  as  Slepian  nodes;  abusing  the  notation 
somewhat,  we  will  retain  this  terminology  while  discussing  higher-dimensional  versions 
of  the  operator  (7),  its  eigenvalues,  and  their  roots. 

We  will  denote  by  An  the  n  x  ra-matrix  defined  by  the  formula 

=  ipj+i(xk),  (9) 

with  Xi,  X21  •  •  • ,  xn  the  roots  of  the  function  ipn  on  the  interval  [—1,1]. 

The  following  theorem  (with  extensions  to  higher  dimensions)  is  the  basis  of  the 
numerical  techniques  we  have  been  developing  recently;  its  proof  is  somewhat  involved 
and  will  be  reported  in  the  paper  currently  in  preparation. 

Theorem  3.2  Suppose  that  for  a  positive  real  a,  xi,X2,“-,xn  are  the  zeroes  on  the 
interval  [—1,1]  of  the  eigenfunction  tpn  of  the  operator  (7),  and  Xn  is  the  eigenvalue 
of  (7)  corresponding  to  rpn.  Then  there  exist  positive  real  numbers  W\,W2 ,•••  ,wn  such 
that 


ll(B")T  ■  Bn  -  /||  <1  An  I, 


where  the  matrix  Bn  is  defined  by  the  formula 


Blk  =  &lk  •  Vk, 


(10) 


(11) 


and  I  denotes  the  unity  matrix.  Furthermore,  for  any  real  0  €  [—2a,  2a], 


Remark  3.2  The  statement  of  Theorem  3.2  can  be  viewed  as  an  extension  of  Gaussian 
quadratures  from  polynomials  to  linear  combinations  of  exponentials,  and  as  a  tool  for  the 
construction  of  interpolation  schemes.  Indeed,  (7)  implies  that  the  nodes  x%,X2 ,  •  •  •  -,xn 
and  corresponding  weights  uq,u>2,  •  •  • ,  wn  integrate  (to  the  precision  |  An  |)  all  functions 
of  the  form  with  (3  €  [—2a,  2a].  Obviously,  the  dimensionality  of  the  space  of 

such  functions  (to  precision  |  An  |)  is  roughly  2  •  n.  In  other  words,  we  have  an  n- 
point  quadrature  with  positive  weights  integrating  accurately  functions  from  a  space  of 
dimensionality  (roughly)  2  •  n. 

Given  a  band-limited  function  /  on  the  interval  [—1, 1]  tabulated  at  the  nodes  aq,  £2,  •  •  • ,  xn, 
an  expansion  of  /  into  a  linear  combination 

/(*)  =  £  7j  *  iM*)  (13) 

j= 0 

is  easily  accomplished  via  the  formula 

7  =  (An)-1(F),  (14) 

with  7  =  (7o,7i>"*,7n-i),  F  =  (f(x0),f{x1),-  ■  •  ,/(xn_i)),  two  vectors  in  Rn.  Since 
the  matrix  Bn  is  nearly  orthogonal  (see  (7)),  the  latter  procedure  is  guaranteed  to  be 
stable  and  provides  a  tool  for  the  interpolation  of  band- limited  functions. 

4  Outline  of  Theory  in  Two  Dimensions 

For  an  arbitrary  real  a  >  0,  we  will  denote  by  Da  the  disk  in  R2  of  radius  a  with 
the  center  at  the  origin.  Given  a  region  Cl  in  R2,  we  will  be  considering  the  operator 
Pa  :  L2(Da)  —y  L2(Cl)  defined  by  the  formula 

PaW){x)  =  [  •  <p(t)dt,  (15) 

JAa 

for  any  ip  €  L2(Da );  the  image  of  L2(Da)  under  Pa  will  be  denoted  by  Xa.  We  will  define 
the  operator  Ba  :  L2(Cl)  — ►  L2(Cl)  by  the  formula 

Ba  =  P*aPa ,  (16) 
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with  P*  denoting  the  adjoint  of  Pa.  A  simple  calculation  shows  that  Ba  is  given  explicitly 
by  the  formula 

(17) 

Obviously,  Ba  is  a  non-negative  definite  symmetric  operator  whose  spectrum  and  eigen¬ 
functions  are  determined  by  the  region  Cl  and  the  parameter  a.  We  will  order  the 
eigenvalues  A,'  of  the  operator  Ba  so  that  Ao  >  Ai  >  A2  >  •  •  •  and  denote  by  1 pi  the 
eigenfunction  corresponding  to  the  eigenvalue  A,-. 

Clearly,  the  operator  Ba  is  compact;  furthermore,  its  spectrum  decays  exponentially. 
An  investigation  of  some  of  the  properties  of  the  operators  Pa,  Ba  can  be  found  in  [4]). 
Given  a  function  f  £  Xa,  it  can  be  represented  in  the  form 

OO 

/(*)  =  £  Si  '  (18) 

j=l 

with  /o,  /i,2  >  •  •  •  a  (unique)  sequence  of  complex  coefficients.  In  a  very  strong  sense,  an 
expansion  of  the  form  (18)  is  the  most  efficient  way  to  represent  a  function  /  €  Xa] 
without  discussing  the  definition  of  an  optimal  representation  in  detail,  we  refer  the 
reader  to  the  classical  paper  [4]  for  a  detailed  discussion  of  such  issues. 

For  each  pair  (i,j)  of  positive  integers  such  that  i  7^  j,  we  will  denote  by  f?tJ  the 
subset  of  Cl  consisting  of  all  points  x  6  Cl  such  that 

Tpi(x)  =  0,  (19) 

M  x)  =  0. 

It  is  easy  to  see  that  for  any  i,j  <  0,  the  set  Rij  consists  of  at  most  a  finite  number  of 
points  in  0;  for  certain  pairs  i,  j,  the  set  Rij  is  empty. 

The  following  observation  is  the  basis  of  the  procedure.  Its  proof  is  involved,  and  will 
be  published  at  a  later  date. 


Observation  4.1  Suppose  that  for  some  integer  k,l  >  0,  the  set  Ritj  consists  ofn  points 
xi,  X2,  •  •  • ,  xn.  Then  there  exist  positive  real  numbers  w\,  w2,  •  •  • ,  wn  such  that 

|  f>r  -  /  e*'^  |<  S(fi)  •  m«*(|Afc|,  |A,|,  |An|),  (20) 

i= i  ,/n 

with  S(Vl)  denoting  the  area  of  Vi),  and  c  an  arbitrary  vector  in  R,2  such  that 

||c||  <2  -  a.  (21) 

In  other  words,  given  two  singular  functions  ipk,  ifi  corresponding  to  eigenvalues  Xk, 
A i,  the  common  zeroes  ofifk,  can  be  viewed  as  the  nodes  of  a  quadrature  formula  on 
f l  for  functions  of  the  form 

e‘-(Ct*),  (22) 

for  all  c  satisfying  the  inequality  (21);  all  weights  of  the  obtained  quadrature  formula  are 
positive.  For  all  such  functions,  the  relative  error  of  the  obtained  quadrature  is  less  than 

max(|Afc|,  |Aj|,  |An|). 

The  following  observation  is  obtained  from  the  preceding  one  via  standard  techniques 
of  approximation  theory. 

Observation  4.2  Suppose  that  under  the  conditions  specified  in  Observation  f.l,  the 
complex  n  x  n-matrix  An  is  defined  by  the  formula 

Alk  =  •  f’j+iixk)-  (23) 

Then 

II (Anr  ■ An-I\\<  S(n)  ■  ma*(|A*|,  |Af|,  |AB|),  (24) 

In  other  words,  the  matrix  connecting  the  coefficients  of  the  (truncated)  expansion  (18), 
with  its  its  values  at  the  points  xx,x2,  -  ••  ,xn  is  almost  orthogonal,  for  sufficiently  large 


Figure  1 

50  nodes  discretizing  20-wavelength  interval  at  3  digits 


5  Examples 

In  this  section,  we  present  three  four  examples  of  receiver  distributions,  two  in  one 
dimension  and  two  in  two  dimensions. 

5.1  Examples  in  One  Dimension 

Both  examples  in  one  dimension  are  of  distributions  of  receivers  on  an  interval  of  a  line. 
In  both  cases,  the  length  of  the  interval  is  20A.  However,  in  the  first  case,  the  relative 
accuracy  of  the  receivers  is  assumed  to  be  1.0F  —  3,  while  in  the  second  case  the  accuracy 
of  receivers  is  assumed  to  be  1.0 E  —  6.  Obviously,  in  the  second  case,  higher  density 
of  receivers  is  needed,  since  the  higher  inherent  sensitivity  of  receivers  will  permit  them 
to  “sense”  the  insufficient  resolution  better.  In  the  first  case  (Figure  1),  the  required 
number  of  receivers  is  50,  or  1.25  times  the  minimum  number  of  nodes  required  by  the 
Nyquist  theorem  (at  fewer  than  2  points  per  wavelength,  an  error  of  100%  is  possible). 
In  the  second  case  (Figure  2),  the  required  number  of  receivers  is  56,  or  1.4  times  the 
Nyquist  minimum. 

5.2  Disks  in  Two  Dimensions 

Here,  we  present  two  examples:  discretization  of  a  disk  of  radius  5  wavelengths  to  3 
digits,  and  discretization  of  the  same  disk  to  6  digits.  As  can  be  seen  in  Figures  3,  4, 


Figure  2 

56  nodes  discretizing  20-wavelength  interval  at  6  digits 


the  5-digit  discretization  requires  677  nodes,  and  the  6-digit  one  requires  796  nodes.  A 
simple  calculation  shows  that  in  the  first  case,  we  expend  about  8.62  nodes  per  square 
wavelength  of  the  region  being  discretized  (about  2.15  times  the  Nyquist  minimum).  In 
the  second  case,  the  corresponding  number  is  10.13  nodes  per  square  wavelength  (about 
2.53  times  the  Nyquist  minimum). 

Remark  5.1  The  algorithms  we  used  to  obtain  results  presented  in  Figures  1-4  are 
based  on  practically  analytical  formulae  for  the  determination  of  the  optimal  nodes; 
these  formulae  are  an  immediate  consequence  of  some  of  the  results  in  [1]-  [4].  The  CPU 
time  requirements  of  the  obtained  numerical  schemes  grow  as  the  square  of  the  length 
of  the  interval  in  one  dimension  (measured  in  the  number  of  wavelengths),  and  as  the 
square  of  the  radios  of  the  disk  in  two  dimensions.  For  all  practical  purposes,  arbitrarily 
large  numbers  of  such  nodes  can  be  constructed  at  a  trivial  cost  in  CPU  time. 

On  the  other  hand,  the  nodes  in  Figures  5-7  are  obtained  by  a  significant  general¬ 
ization  of  the  results  of  [4];  the  resulting  algorithms  are  extremely  inefficient  (the  cost 
grows  as  the  6-th  power  of  the  size  of  the  triangle!).  The  situation  is  salvaged  to  a  large 
extent  by  the  fact  that  each  configuration  has  to  be  constructed  only  once,  so  that  large 
expenditure  of  the  CPU  time  is  often  affordable.  Still,  triangles  in  Figures  5-7  are  much 
smaller  than  the  disk  in  Figures  3,4. 


Figure  3 

677  nodes  discretizing  to  three  digits  a  disk  of  radius  5  wavelengths 


Figure  4 

796  nodes  discretizing  to  six  digits  a  disk  of  radius  5  wavelengths 
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5.3  Triangles  in  Two  Dimensions 

Here,  we  present  three  examples:  discretizations  of  an  equilateral  triangle  with  the  side 
of  length  1A  to  3  and  5  digits,  and  a  discretization  of  an  equilateral  triangle  with  the 
side  of  length  2A  to  3  digits.  As  can  be  seen  in  Figures  5,  6,  the  3-digit  discretization 
requires  22  nodes,  and  the  6-digit  one  requires  40  nodes.  The  first  case  corresponds  to 
about  51  nodes  per  square  wavelength,  and  the  second  corresponds  to  about  92  nodes 
per  square  wavelength.  A  3-digit  discretization  of  a  triangle  with  the  side  of  length  2A 
is  depicted  in  In  Figure  7.  It  requires  40  nodes,  which  ttranslates  into  about  23  nodes 
per  square  wavelength. 

Remark  5.2  At  first  glance,  it  seems  strange  that  the  discretization  in  Figure  7  requires 
roughly  twice  as  many  nodes  as  that  in  Figure  6,  while  the  area  of  the  triangle  is  4  times 
greater.  However,  this  behavior  is  not  at  all  anomalous,  since  for  a  sufficiently  large 
region,  the  required  number  of  nodes  should  converge  to  4  points  per  square  wavelength, 
as  required  by  the  Nyquist  theorem.  Generally,  the  required  density  of  nodes  increases 
with  the  decrease  in  the  size  of  the  region  being  discretized. 

6  Application  to  the  Measurements  and  Generation 
of  Band-Limited  Data 

Due  to  the  source-receiver  duality  in  the  antenna  theory,  we  will  be  only  discussing 
the  applications  of  the  techniques  in  the  design  of  receiving  antennae.  The  situation  is 
identical  in  the  design  of  radiating  ones. 

Assuming  that  our  receivers  are  distributed  on  a  specified  region  Cl  in  R2,  and  their 
sensitivity  e  is  specified,  we  would  like  to  find  the  distribution  of  the  receivers  that  would 
obtain  the  accuracy  e  in  the  measurement  of  any  function  of  the  form  (22),  with  ||c||  <  a. 

Our  procedure  calls  for  constructing  (numerically)  the  eigenvalues  and  eigenvectors  of 
the  operator  Ba ,  and  examining  the  common  zeroes  of  several  pairs  of  eigenvectors  of  Ba 


Figure  5 

22  nodes  discretizing  to  3  digits  an  equilateral  triangle  with  the  side  1A 


Figure  6 

40  nodes  discretizing  to  5  digits  an  equilateral  triangle  with  the  side  1A 


Figure  7 

40  nodes  discretizing  to  3  digits  an  equilateral  triangle  with  the  side  2A 

whose  corresponding  eigenvalues  are  close  to  e.  Supposing  that  two  such  eigenfunctions 
tjV’i  have  n  zeroes  in  common,  and  that  |A*|  <  e,  |A/|  <  e,  |An|  <  e,  the  obtained 
nodes  provide  a  satisfactory  set  of  locations  for  the  antennae.  Choosing  among  such 
satisfactory  sets  the  one  for  which  n  is  minimized,  we  obtain  the  optimum  distribution 
of  antennae  for  the  accuracy  e  and  region  Cl. 

Remark  6.1  Given  a  region  in  the  plane  (or  on  a  surface  in  R3)  where  the  phased 
array  antenna  is  to  be  located,  the  analytical  apparatus  we  have  constructed  provides 
(at  some  computational  cost)  optimal  locations  for  the  receivers  or  transducers.  Clearly, 
this  does  not  address  the  many  engineering  issues  that  will  arize  in  connection  with  such 
configurations.  At  the  present  time,  we  are  investigating  these  issues,  while  at  the  same 
time  preparing  a  patent  application  for  the  procedure. 


7  Application  to  the  Discretization  of  Surfaces  of 
Scatterers 

When  the  technique  is  to  be  applied  to  the  design  of  discretizations  of  surfaces  of  scatter¬ 
ers  in  the  numerical  modeling  of  electromagnetic  and  acoustic  phenomena,  we  start  with 
a  triangulation  of  the  surface.  We  will  assume  that  the  discretization  provided  to  us  sat¬ 
isfies  the  standard  condition  that  the  triangles  on  the  surface  have  limited  aspect  ratios. 
It  is  usual  to  parametrize  the  triangles  in  a  triangulation  by  the  standard  equilateral 
triangle  in  R2,  and  we  assume  that  our  triangles  are  so  parametrized. 

Given  a  user-specified  accuracy  and  the  size  of  the  triangle  on  the  surface  (we  mea¬ 
sure  the  size  in  the  number  of  wavelengths  of  the  incoming  radiation),  we  construct  an 
optimum  discretization  of  the  standard  triangle  in  the  plane  that  guarantees  the  required 
precision  at  a  frequency  somewhat  higher  than  that  at  which  the  problem  is  being  solved. 
The  actual  amount  of  oversampling  is  determined  by  the  smoothness  of  the  surface  in 
the  region  being  discretized;  most  of  the  time,  it  will  be  determined  predominantly  by 
the  quality  of  the  CAD  package  that  produced  the  surface  being  modeled.  Now,  the 
images  on  the  surface  of  the  scatterer  of  the  nodes  on  the  standard  triangle  are  the 
nodes  of  the  discretization  of  the  surface.  The  associated  quadrature  and  interpolation 
formulae  on  the  standard  triangle  become  (with  appropriate  scaling)  the  quadrature  and 
interpolation  formulae  on  the  surface. 

Remark  7.1  It  should  be  observed  that  while  the  theory  of  such  discretizations  is  quite 
straightforward  (given  the  discretization  of  the  standard  triangle),  we  have  only  tested 
the  resulting  numerical  procedure  in  two  dimensions.  Obviously,  several  numerical  issues 
have  to  be  resolved  before  the  approach  becomes  an  effective  tool  in  the  numerical 
solution  of  scattering  problems  in  three  dimensions. 
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A  Procedure  for  the  Design  of  Apparata  for  the  Measurement 

of  Band-Limited  Signals 


Abstract 

Whenever  physical  signals  are  measured  or  generated,  the  locations  of  receivers  or 
transducers  have  to  be  selected.  Most  of  the  time,  this  appears  to  be  done  on  an  ad  hoc 
basis.  For  example,  when  a  string  of  geophones  is  used  in  the  measurements  of  seismic 
data  in  oil  exploration,  the  receivers  are  located  at  equispaced  points  on  an  interval. 
When  phased  array  antennae  are  constructed,  their  shapes  are  determined  by  certain 
aperture  considerations;  round  and  rectangular  shapes  are  common.  When  antenna 
beams  are  steered  electronically,  it  is  done  by  changing  the  phases  and  the  amplitudes 
of  of  the  transducers.  Again,  these  transducers  are  located  in  a  region  of  predetermined 
geometry,  and  their  actual  locations  within  that  geometry  are  chosen  via  some  heuristic 
procedure. 

In  all  these  (and  many  other)  cases,  the  signals  being  received  or  generated  are  band- 
limited.  Optimal  representation  of  such  signals  has  been  studied  in  detail  in  [l]-[4]  more 
than  30  years  ago.  We  use  results  from  [l]-[4]  (with  some  extensions)  to  construct  optimal 
nodes  for  the  measurement  and  generation  of  band-limited  signals.  In  this  disclosure, 
we  describe  the  procedure  based  on  these  techniques  for  the  design  of  such  receiver  (and 
transducer)  configurations  in  a  variety  of  environments. 

1  Introduction 

When  measurements  are  performed,  it  often  happens  that  the  signal  to  be  measured  is 
well  approximated  by  linear  combinations  of  oscillatory  exponentials,  i.e.  functions  of 


(2) 


in  one  dimension,  of  the  form 

•  ei-{x^x+^y) 
i= i 

in  two  dimensions,  and  of  the  form 

n 

•  ei'{'X]'x+y'3'yJrl'3'z^  (3) 

i=i 

in  three  dimensions.  In  most  cases,  the  signal  is  band-limited,  i.e.  there  exist  such  real 
positive  a  that  all  1  <  j  <  n, 

I  A j  |<  a  (4) 

in  one  dimension, 

Aj  4-  nj  <  a2  (5) 

in  two  dimensions,  and 

A|  +  /^ |  +  v2  ^  a2,  (6) 

in  three  dimensions. 

As  is  well-known,  most  measurements  of  electromagnetic  and  acoustic  data  (espe¬ 
cially  at  reasonably  high  frequencies)  are  of  this  form.  Examples  of  such  situations 
include  geophone  and  hydrophone  strings  in  geophysics,  phased  array  antennae  in  radar 
systems,  multiple  transceivers  in  ultrasound  imaging,  and  a  number  of  other  applications 
in  astrophysics,  medical  imaging,  non-destructive  testing,  etc. 

In  this  disclosure,  we  describe  a  procedure  for  determining  the  optimal  distribution  of 
sources  and  receivers  that  maximizes  accuracy  and  resolution  in  measuring  band-limited 
data  given  a  fixed  number  of  receivers.  Alternatively,  the  procedure  can  be  used  to 
determine  the  optimal  distribution  of  receivers  that  will  minimize  their  number  given 
specified  accuracy  and  resolution.  We  illustrate  the  procedure  with  several  examples 
of  optimal  distributions  of  receivers  on  an  interval  (one-dimensional  case),  on  the  disk 
(two-dimensional  case),  and  on  the  triangle  (two-dimensional  case). 
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Remark  1.1  Our  procedure  permits  the  design  of  optimal  receiver  distributions  in  ar¬ 
bitrary  regions  in  one,  two,  and  three  dimensions.  However,  when  the  region  in  question 
is  an  interval,  a  disk,  or  a  sphere,  there  exists  an  explicit  formula  for  the  locations  of 
the  receivers,  in  terms  of  zeroes  of  the  Generalized  Prolate  Spheroidal  Wave  Functions. 
Somewhat  more  complicated,  but  still  nearly  analytical,  schemes  exist  when  the  region 
is  the  triangle  (in  two  dimensions),  a  tetrahedron  (in  three  dimensions),  and  in  several 
other  cases.  In  order  to  avoid  describing  the  numerical  apparatus  required  in  the  general 
case,  we  restrict  our  examples  to  the  cases  where  the  problem  can  be  solved  explicitly 
(as  opposed  to  an  involved  numerical  procedure). 

The  analytical  apparatus  we  use  is  closely  related  both  to  classical  Gaussian  quadra¬ 
tures  and  to  the  analysis  of  band- limited  functions  initiated  in  the  classical  works  [l]-[4]. 
It  results  in  optimal  (in  a  very  strong  sense)  choice  for  both  source  and  receiver  loca¬ 
tions.  While  the  construction  is  quite  straightforward  in  one  dimension  (in  this  case,  it 
is  a  remarkably  simple  corollary  to  [1],  [2]),  in  higher  dimensions  it  is  complicated  by 
the  fact  that  regions  can  be  of  various  shapes.  In  practical  measurements,  the  cases  of 
the  string,  of  the  disk,  of  the  square,  and  of  the  triangle  appear  to  be  most  frequently 
encountered. 

2  Outline  of  Theory  In  One  Dimension 

Given  a  real  a  >  0,  we  will  be  considering  the  operator  Pa  :  L2[— 1, 1]  — *  L2[— 1, 1]  defined 
by  the  formula 

(7) 

for  any  tp  €  L2[— 1,1];  the  image  of  L2  under  Pa  will  be  denoted  by  Xa.  As  observed 
in  [1],  all  eigenvalues  A0,  Ai,  A2,  •  •  •  of  the  operator  Pa  are  distinct,  and  the  corresponding 
eigenvectors  are  the  appropriately  chosen  Prolate  Spheroidal  Wave  functions. 

The  following  theorem  is  a  combination  of  several  Lemmas  from  [1]. 


Theorem  2.1  For  any  positive  real  a,  the  eigenfunctions  ip o,ipi,---,  of  the  operator 
Pa  constitute  an  orthonormal  basis  in  Xa.  The  function  ipj  has  exactly  i  zeroes  on  the 
interval  [—1, 1],  and  for  any  positive  integer  n,  the  functions  ipo ,  ipi,  ■  •  •  ,ipn  constitute  a 
Chebychev  system  on  the  interval  [—1,1]. 

Assuming  that  A0,Ai,A2,---  are  ordered  so  that  |  Ao  |>|  Ai  |>|  A2  |>  •••,  the  even- 
numbered  eigenvalues  are  purely  real,  and  the  odd-numbered  ones  are  purely  imaginary. 
The  even-numbered  eigenfunctions  are  even,  the  odd-numbered  ones  are  odd,  and  the 
j  —  th  eigenfunction  has  exactly  j  roots  on  the  interval  [—1, 1]. 

Finally,  for  all  sufficiently  large  i, 

I  A;  |<  e-i.  (8) 

Remark  2.1  The  operator  (7)  and  some  of  its  generalizations  have  been  studied  in  [1]- 
[5];  later,  a  systematic  investigation  of  generalizations  of  (7)  was  undertaken  in  [7]- [9]. 
The  investigations  reported  in  [l]-[5]  appear  to  have  been  motivated  almost  entirely  by 
considerations  of  Electrical  Engineering;  the  work  [7]- [9]  seems  to  have  been  motivated 
predominantly  by  considerations  of  pure  mathematics.  Remarkably,  implications  of  this 
mathematical  apparatus  in  numerical  analysis  appear  to  have  been  overlooked. 

In  the  interests  of  brevity,  we  will  be  referring  to  the  eigenfunctions  ipn  of  the  oper¬ 
ator  (7)  as  Slepian  functions,  and  to  their  roots  as  Slepian  nodes;  abusing  the  notation 
somewhat,  we  will  retain  this  terminology  while  discussing  higher- dimensional  versions 
of  the  operator  (7),  its  eigenvalues,  and  their  roots. 

We  will  denote  by  An  the  n  x  n-matrix  defined  by  the  formula 

Alk  =  (9) 

with  ®i,  X2,  •  •  • ,  xn  the  roots  of  the  function  ipn  on  the  interval  [—1,1]. 

The  following  theorem  (with  extensions  to  higher  dimensions)  is  the  basis  of  the 
numerical  techniques  we  have  been  developing  recently;  its  proof  is  somewhat  involved 
and  will  be  reported  in  the  paper  currently  in  preparation. 


Theorem  2.2  Suppose  that  for  a  positive  real  a,  ®i,  X2>  • '  ‘  >  xn  are  the  zeroes  on  the 
interval  [—1,1]  of  the  eigenfunction  f>n  of  the  operator  (7),  and  An  is  the  eigenvalue 
of  (7)  corresponding  to  f>n.  Then  there  exist  positive  real  numbers  W\,  iw2,  •  •  • ,  wn  such 
that 


\\(Bn)T  ■  B"  —  /||  <|  A„  |, 


where  the  matrix  Bn  is  defined  by  the  formula 


(10) 


Blk  =  Alk 


Wk, 


(11) 


and  I  denotes  the  unity  matrix.  Furthermore,  for  any  real  j3  €  [—2a,  2a], 

lX>j-  -  I'  e'^dx  |<|  \n  |,  (12) 

i=i  J~1 

Remark  2.2  The  statement  of  Theorem  2.2  can  be  viewed  as  an  extension  of  Gaussian 
quadratures  from  polynomials  to  linear  combinations  of  exponentials,  and  as  a  tool  for  the 
construction  of  interpolation  schemes.  Indeed,  (7)  implies  that  the  nodes  xi,x2,  •  • • , xn 
and  corresponding  weights  u>i,u;2,  •  •  •  ,wn  integrate  (to  the  precision  j  A„  |)  all  functions 
of  the  form  el'®'Xi ,  with  (3  €  [—2a,  2a].  Obviously,  the  dimensionality  of  the  space  of 
such  functions  (to  precision  |  An  |)  is  roughly  2  •  n.  In  other  words,  we  have  an  n- 
point  quadrature  with  positive  weights  integrating  accurately  functions  from  a  space  of 
dimensionality  (roughly)  2  •  n. 

Given  a  band-limited  function  /  on  the  interval  [—1,1]  tabulated  at  the  nodes  xi,  x2,  •  •  • , 
an  expansion  of  /  into  a  linear  combination 

fix)  =  £  7j  •  Mx)  (13) 

3=0 

is  easily  accomplished  via  the  formula 

7  =  (xr‘(n  (w) 

with  7  =  (7o57i)*"  )7n-i)i  F  —  ,  f{Xn-i )),  two  vectors  in  It".  Since 

the  matrix  Bn  is  nearly  orthogonal  (see  (7)  ),  the  latter  procedure  is  guaranteed  to  be 
stable  and  provides  a  tool  for  the  interpolation  of  band-limited  functions. 


3  Outline  of  Theory  in  Two  Dimensions 


For  an  arbitrary  real  a  >  0,  we  will  denote  by  Da  the  disk  in  R2  of  radius  a  with 
the  center  at  the  origin.  Given  a  region  0  in  R,2,  we  will  be  considering  the  operator 
Pa  :  L2(Da )  — ►  L2(Vt)  defined  by  the  formula 

Pa(<p)(x)  =  f  (15) 

J Aa 

for  any  tp  €  L2(Da)]  the  image  of  L2(Da )  under  Pa  will  be  denoted  by  Xa.  We  will  define 
the  operator  Ba  :  L2{Vt)  — ►  L2( fi)  by  the  formula 

Ba  =  Pa*P0,  (16) 

with  P*  denoting  the  adjoint  of  Pa.  A  simple  calculation  shows  that  Ba  is  given  explicitly 
by  the  formula 

Ba(v)(x)  =  Jq  l(||JLt||  ~dt  (1?) 

Obviously,  Ba  is  a  non-negative  definite  symmetric  operator  whose  spectrum  and  eigen¬ 
functions  are  determined  by  the  region  f l  and  the  parameter  a.  We  will  order  the 
eigenvalues  A,-  of  the  operator  Ba  so  that  Ao  >  Ai  >  A2  >  •  •  •  and  denote  by  the 
eigenfunction  corresponding  to  the  eigenvalue  A,\ 

Clearly,  the  operator  Ba  is  compact;  furthermore,  its  spectrum  decays  exponentially. 
An  investigation  of  some  of  the  properties  of  the  operators  Pa ,  Ba  can  be  found  in  [4]). 
Given  a  function  f  (E  Xa,  it  can- be  represented  in  the  form 

/(*)  =  £  fj  •  &(*)»  (18) 

i=l 

with  /o,  /i,2  ,  •  •  •  a  (unique)  sequence  of  complex  coefficients.  In  a  very  strong  sense,  an 
expansion  of  the  form  (18)  is  the  most  efficient  way  to  represent  a  function  f  €  X 0; 
without  discussing  the  definition  of  an  optimal  representation  in  detail,  we  refer  the 
reader  to  the  classical  paper  [4]  for  a  detailed  discussion  of  such  issues. 


For  each  pair  (i,j)  of  positive  integers  such  that  i  ^  j,  we  will  denote  by  Rij  the 
subset  of  Cl  consisting  of  all  points  x  6  Cl  such  that 

Mx)  =  (19) 

i>j{x)  =  o. 

It  is  easy  to  see  that  for  any  i,j  <  0,  the  set  Rij  consists  of  at  most  a  finite  number  of 
points  in  Cl ;  for  certain  pairs  i,j,  the  set  Rij  is  empty. 

The  following  observation  is  the  basis  of  the  procedure  being  disclosed.  Its  proof  is 
involved,  and  will  be  published  at  a  later  date. 


Observation  3.1  Suppose  that  for  some  integer  k,l  >  0,  the  set  Rij  consists  of  n  points 
X\, #2, '  ■ "  j 3?n •  Then  there  exist  positive  real  numbers  wi,W2 ,  •  •  • ,  wn  such  that 

IX>r  |<  S(fi)  •  mo*(|A»|,  | Ai|,  |A,|),  (20) 

fr.  Ja 

with  S{Cl)  denoting  the  area  of  Cl),  and  c  an  arbitrary  vector  in  R,2  such  that 

||c||  <2  -  a.  (21) 


In  other  words ,  given  two  singular  functions  xfk,  tfi  corresponding  to  eigenvalues  A*, 
Af,  the  common  zeroes  ofifk,  can  be  viewed  as  the  nodes  of  a  quadrature  formula  on 
Cl  for  functions  of  the  form 


(22) 


for  all  c  satisfying  the  inequality  (21);  all  weights  of  the  obtained  quadrature  formula  are 
positive.  For  all  such  functions,  the  relative  error  of  the  obtained  quadrature  is  less  than 
max(|Afc|,|A/|,|An|). 


The  following  observation  is  obtained  from  the  preceding  one  via  standard  techniques 
of  approximation  theory. 


Observation  3.2  Suppose  that  under  the  conditions  specified  in  Observation  3.1,  the 
complex  n  x  n-matrix  An  is  defined  by  the  formula 

Aj,k  =  Vwk  •  ipj+i(xk).  (23) 

Then 

|| (Any  •  An  -  /||  <  S(Q)  ■  max( |A*|, |A,|,  |An|),  (24) 

In  other  words,  the  matrix  connecting  the  coefficients  of  the  (truncated)  expansion  (18), 
with  its  its  values  at  the  points  X\,X2,  •  •  *  is  almost  orthogonal,  for  sufficiently  large 
k ,  l,n. 

4  Examples 

In  this  section,  we  present  three  four  examples  of  receiver  distributions,  two  in  one 
deimension  and  two  in  two  dimensions. 

4.1  Examples  in  One  Dimension 

Both  examples  in  one  dimension  are  of  distributions  of  receivers  on  an  interval  of  a  line. 
In  both  cases,  the  length  of  the  interval  is  20A.  However,  in  the  first  case,  the  relative 
accuracy  of  the  receivers  is  assumed  to  be  1.0E  —  3,  while  in  the  second  case  the  accuracy 
of  receivers  is  assumed  to  be  l.OE  —  6.  Obviously,  in  the  second  case,  higher  density 
of  receivers  is  needed,  since  the  higher  inherent  sensitivity  of  receivers  will  permit  them 
to  “sense”  the  insufficient  resolution  better.  In  the  first  case  (Figure  1),  the  required 
number  of  receivers  is  50,  or  1.25  times  the  minimum  number  of  nodes  required  by  the 
Nyquist  theorem  (at  fewer  than  2  points  per  wavelength,  an  error  of  100%  is  possible). 
In  the  second  case  (Figure  2),  the  required  number  of  receivers  is  56,  or  1.4  times  the 
Nyquist  minimum. 


Figure  1 

50  nodes  discretizing  20-wavelength  interval  at  3  digits 


Figure  2 

56  nodes  discretizing  20- wavelength  interval  at  6  digits 


4.2  Disks  in  Two  Dimensions 


Here,  we  present  two  examples:  discretization  of  a  disk  of  radius  5  wavelengths  to  3 
digits,  and  discretization  of  the  same  disk  to  6  digits.  As  can  be  seen  in  Figures  3,  4, 
the  5-digit  discretization  requires  677  nodes,  and  the  6-digit  one  requires  796  nodes.  A 
simple  calculation  shows  that  in  the  first  case,  we  expend  about  8.62  nodes  per  square 
wavelength  of  the  region  being  discretized  (about  2.15  times  the  Nyquist  minimum).  In 
the  second  case,  the  corresponding  number  is  10.13  nodes  per  square  wavelength  (about 
2.53  times  the  Nyquist  minimum). 

Remark  4.1  The  algorithms  we  used  to  obtain  results  presented  in  Figures  1-4  are 
based  on  practically  analytical  formulae  for  the  determination  of  the  optimal  nodes; 
these  formulae  are  an  immediate  consequence  of  some  of  the  results  in  [1]-  [4].  The  CPU 
time  requirements  of  the  obtained  numerical  schemes  grow  as  the  square  of  the  length 
of  the  interval  in  one  dimension  (measured  in  the  number  of  wavelengths),  and  as  the 
square  of  the  radios  of  the  disk  in  two  dimensions.  For  all  practical  purposes,  arbitrarily 
large  number  oif  such  nodes  can  be  construcrted  at  a  trivial  cost  in  CPU  time. 

On  the  other  hand,  the  nodes  in  Figures  5-7  are  obtained  by  a  significant  general¬ 
ization  of  the  results  of  [4];  the  resulting  algorithms  are  extremely  inefficient  (the  cost 
grows  as  the  6-th  power  of  the  size  of  the  triangle!).  The  situation  is  salvaged  to  a  large 
extent  by  the  fact  that  each  configuration  has  to  be  constructed  only  once,  so  that  large 
expenditure  of  the  CPU  time  is  often  affordable.  Still,  triangles  in  Figures  5-7  are  much 
smaller  than  the  disk  in  Figures  3,4. 

4.3  Triangles  in  Two  Dimensions 

Here,  we  present  three  examples:  discretizations  of  an  equilateral  triangle  with  the  side 
of  length  1A  to  3  and  5  digits,  and  a  discretization  of  an  equilateral  triangle  with  the 
side  of  length  2 A  to  3  digits.  As  can  be  seen  in  Figures  5,  6,  the  3-digit  discretization 


Figure  3 

677  nodes  discretizing  to  three  digits  a  disk  of  radius  5  wavelengths 


Figure  4 

796  nodes  discretizing  to  three  digits  a  disk  of  radius  5  wavelengths 
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Figure  5 

22  nodes  discretizing  to  3  digits  an  equilateral  triangle  with  the  side  1A 

requires  22  nodes,  and  the  6-digit  one  requires  40  nodes.  The  first  case  corresponds  to 
about  51  nodes  per  squaxe  wavelength,  and  the  second  corresponds  to  about  92  nodes 
per  square  wavelength.  A  3-digit  discretization  of  a  triangle  with  the  side  of  length  2A  is 
depicted  in  In  Figure  7.  It  requires  40  nodes,  which  translates  into  about  23  nodes  per 
square  wavelength. 

Remark  4.2  At  first  glance,  it  seems  strange  that  the  discretization  in  Figure  7  requires 
roughly  twice  as  many  nodes  as  that  in  Figure  6,  while  the  area  of  the  triangle  is  4  times 
greater.  However,  this  behavior  is  not  at  all  anomalous,  since  for  a  sufficiently  large 
region,  the  required  number  of  nodes  should  converge  to  4  points  per  square  wavelength, 
as  required  by  the  Nyquist  theorem.  Generally,  the  required  density  of  nodes  increases 
with  the  decrease  in  the  size  of  the  region  being  discretized. 

5  Applications  to  the  Measurements  and  Generation 
of  Band-Limited  Data 

Due  to  the  source-receiver  duality  in  the  antenna  theory,  we  will  be  only  discussing 
the  applications  of  the  techniques  in  the  design  of  receiving  antennae.  The  situation  is 


Figure  6 

40  nodes  discretizing  to  5  digits  an  equilateral  triangle  with  the  side  1A 


Figure  7 

40  nodes  discretizing  to  3  digits  an  equilateral  triangle  with  the  side  2A 


identical  in  the  design  of  radiating  ones. 

Assuming  that  our  receivers  are  distributed  on  a  specified  region  0  in  K,2,  and  their 
sensitivity  e  is  specified,  we  would  like  to  find  the  distribution  of  the  receivers  that  would 
obtain  the  accuracy  e  in  the  measurement  of  any  function  of  the  form  (22),  with  ||c||  <  a. 

Our  procedure  calls  for  constructing  (numerically)  the  eigenvalues  and  eigenvectors  of 
the  operator  Ba,  and  examining  the  common  zeroes  of  several  pairs  of  eigenvectors  of  Ba 
whose  corresponding  eigenvalues  are  close  to  e.  Supposing  that  two  such  eigenfunctions 
have  n  zeroes  in  common,  and  that  | A* |  <  e,  |A;|  <  e,  |An|  <  e,  the  obtained 
nodes  provide  a  satisfactory  set  of  locations  for  the  antennae.  Choosing  among  such 
satisfactory  sets  the  one  for  which  n  is  minimized,  we  obtain  the  optimum  distribution 
of  antennae  for  the  accuracy  e  and  region  fi. 
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